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Abstract. Nonlinear Schrodinger Equations (NLS) of the Hartree type occur 
in the modeling of quantum semiconductor devices. Their " semiclassical" limit 
of vanishing (scaled) Planck constant is both a mathematical challenge and 
practically relevant when coupling quantum models to classical models. With 
the aim of describing the semi-classical limit of the 3D Schrodingcr-Poisson 
system with an additional harmonic potential, we study some semi-classical 
limits of the Hartree equation with harmonic potential in space dimension 
n > 2. The harmonic potential is confining, and causes focusing periodically 
in time. We prove asymptotics in several cases, showing different possible non- 
linear phenomena according to the interplay of the size of the initial data and 
the power of the Hartree potential. In the case of the 3D Schrodinger— Poisson 
system with harmonic potential, we can only give a formal computation since 
the need of modified scattering operators for this long range scattering case 
goes beyond current theory. 

We also deal with the case of an additional "local" nonlinearity given by a 
power of the local density - a model that is relevant when incorporating the 
Pauli principle in the simplest model given by the " Schrodinger-Poisson-Xa 
equation" . Further we discuss the connection of our WKB based analysis to 
the Wigner function approach to semiclassical limits. 



1. Introduction 

Nonlinear Schrodinger Equations (NLS) are important both for many differ- 
ent applications as well as a source of rich mathematical theory with several hard 
challenges still open. The NLS in the most common meaning contains a "local" 
nonlinearity given by a power of the local density, in particular the (de)focusing 
"cubic" NLS which arises e.g. in nonlinear optics or for Bose Einstein condensates. 
In 1-d this NLS is an integrable system and the "semi-classical limit" ("high wave 
number limit") can be performed by methods of inverse scattering (see e.g. (201 an d 
[22] for results on the defocusing and focusing case). A class of NLS with a "non- 
local" nonlinearity that we call "Hartree type" occur in the modeling of quantum 
semiconductor devices. Their "semi-classical" limit of vanishing (scaled) Planck 
constant is both a mathematical challenge and practically relevant when coupling 
quantum models to classical models. 

Incorporating the Pauli principle for fermions in the simplest possible model yields 
the case of a Hartree equation with an additional "local" nonlinearity given by a 
power of the local density, the " Schrodinger-Poisson-Xa equation" (see [5K]L 

In this paper we deal with the "semi-classical limit" of nonlinear Schrodinger 
equations of Hartree type, with a harmonic potential and a "weak" nonlinearity 
which is a convolution of the density with a more or less singular potential. 
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In three space dimensions, for the case where we convolute with the Newto- 
nian potential l/|x|, the Hartree equation is the Schrodinger-Poisson system with 
harmonic potential : 

ied t u e + ^e 2 An e = ^u £ + V{x)u £ , 
AV=\u £ \ 2 , 
u ft=o = u ) 
with x € M 3 . 

This equation arises typically if we consider the quantum mechanical time evolu- 
tion of electrons in the mean field approximation of the many body effects, modeled 
by the Poisson equation, with a confinement modeled by the quadratic potential of 
the harmonic oscillator. 

The limit e — > in such a quantum model corresponds to a "classical limit" of 
vanishing Planck constant h = e — > 0. We adopt the terminology "semi-classical 
limit" for what should properly be called "classical limit" (see the discussion in 
the term "semi-classical" being actually more appropriate for the situation of 
the homogenization limit from a Schrodinger equation with periodic potential (see 
e.g. 0). 

The problem of the mathematically rigorous "classical limit" of the Schrodinger- 
Poisson system is highly nontrivial. First results of weak limits e — ■> to the Vlasov- 
Poisson system where given in |2,'ij and |24| using Wigner transform techniques for 
the "mixed state case" , where additional strong assumptions on the initial data 
can be imposed (which are necessary to guarantee a uniform L 2 bound on the 
Wigner function). In [31] this assumption could be removed for the 1-d case and 
the classical limit for the "pure state" case could be performed, where the notorious 
problem of non-uniqueness of the Vlasov-Poisson system with measure valued initial 
data reappears. For an overview of this kind of "semi-classical limits" of Hartree 
equations see |26) . For an introduction to Wigner transforms and their comparison 
to WKB methods for the linear case see and |29j . 

Up to a constant, (|l.l|l is equivalent to the Hartree equation 

(1.2) ?£ c^ + ie 2 Au^^ + (M-^K| 2 )u e ; 4=o = u o- 

We restrict our attention to small data cases with = e a ^ 2 f, where / is indepen- 
dent of e and a > 1. 

Notice that we can allow for more general data with initial plane oscillations, 

(1.3) uf t=0 = s a ' 2 f(xy^ for &€M 3 , 
since the change of variables given in |H] 

(1.4) u E (t,x) =u e (t,x~£ smt)e< x -^ sint )-t° cost/e , 

yields the solution of 1(1. 2(1 . This change of variable could also be used in Equa- 
tion 1(1.6(1 below and hence our results also hold for the more general e-dependent 
class of data (|1.3fl . 

Note that "small data" can be equivalently written as "small nonlinearity" , since 
with the change of the unknown u e = e _Q / 2 u e , p. 2(1 becomes 

(1.5) ied t u s + ^e 2 Au e = ^-u s + e a (\x\- 1 *\u e \ 2 )u s ; uf t=0 = /. 
We will consider the more general "semi-classical Hartree equation" 

(1.6) zed t u £ + ^e 2 Au^^u' + e a (\x\~^\u'\ 2 )u' ; uf t=0 = /, 
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with 7 > 0, a > 1 and x € M. n , where the space dimension n > 2 may be different 
from 3. 

The first point to notice is that in the linear case, the harmonic potential causes 
focusing at the origin (resp. at (— l) fe £o m the case (II .4(1 ) at times t = ir/2 + kir, 
for any k G N. The solution itf rce of the linear equation 

1 \x\ 2 

(1.7) i£d t uf Iee + -£ 2 A?4 ec = —u% ee ; < cc |t=o = / » 

is initially of size 0(1). At time t — tt/2, the solution focuses at the origin and is of 
order 0(e - "/ 2 ); it is of order 0(1) for t = 7r, and so on (for a more precise analysis, 
see |H|). This phenomenon is easy to read from Mehler's formula (see e.g. |l()M18j ): 
for < t < ir, we have 

Essentially, one can apply a stationary phase formula for t G]0, 7r/2[U]7r/2, tt[ (uf rcc 
is 0(1)), while it is not possible at t — n/2 (itf ree is 0(e~™/ 2 )). Following the same 
approach as in [3], we get the following distinctions: 





a > 7 


a — 7 


a > 1 


Linear WKB, 


Linear WKB, 




linear focus 


nonlinear focus 


a = 1 


Nonlinear WKB, 


Nonlinear WKB, 




linear focus 


nonlinear focus 



The expression "linear WKB" means that the nonlinear Hartree interaction term 
is negligible away from the focus (when the WKB approximation is valid); "linear 
focus" means that the nonlinearity is negligible near the focus; the WKB regime 
(resp. the focus) is "nonlinear" when the Hartree term has a leading order influence 
away from (resp. in the neighborhood of) the focus, in the limit s — > 0. This 
terminology follows |19j. 

We did not obtain a rigorous description of the case a = 7 = 1, which corre- 
sponds to the Schrodinger-Poisson system (|l.lfl when n = 3. This problem seems 
out of reach for the methods currently available in this held. On the other hand, 
we study rigorously the three other cases in an exhaustive way: 

In Section we prove that the Hartree term has no influence at leading order 
when a > 7 = 1. Back to Q1.2JI . this shows that initial data of size e Q / 2 with 
a > 1 yield a linearizable solution. The expected critical size is y/e; this heuristic 
is reinforced by the next three sections. 

In Section we study the case a = 1 > 7. We prove that the nonlinear term 
must be taken into account to describe the solution u £ . It is so through a slowly 
oscillating phase term. On the other hand, no nonlinear effect occurs at leading 
order near the focus. 

In Section |5l we show that when a = 7 > 1, nonlinear effects occur at leading 
order at the focuses, while they are negligible elsewhere. This phenomenon is 
the same as in for the nonlinear Schrodinger equation; each focus crossing is 
described in terms of the scattering operator associated to the Hartree equation 

(1.9) idtili + ±Ail>= * \ip\ 2 ) . 

In Sectional we perform a formal computation suggested by the results of Sec- 
tions 0| and |SJ This can be seen as a further evidence that nonlinear effects are 
always relevant in the case a = 7 = 1, along with a precise idea of the nature of 
these nonlinear effects, which we expect to be true. We add a brief discussion of 
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the case of an additional local nonlinearity in the equation and some remarks on 
the Wigner measures in view of the ill-posedness results of [H] • 

This program is very similar to the one achieved in 0. We want to underscore 
at least two important differences. First, we have to adapt the notion of oscillatory 
integral to incorporate the presence of the harmonic potential (see Section 14.10 . 
Second, the power-like nonlinearity treated in [3] is replaced by a Hartree-type 
nonlinearity. This yields different and less technical proofs (we do not use Strichartz 
estimates in Sections and 0J, and makes a more complete description of the 
above table possible; the case "nonlinear WKB, linear focus" was treated very 
partially in [3], due to the lack of regularity of the map z i— ► \z\ 2a z for small a > 0. 
This technical difficulty does not occur in the present case, and the main result of 
Section EH (Proposition |4~TJ> is proved with no restriction. 

The content of this article is as explained above, plus a paragraph dedicated to a 
quick review of the facts we will need about the Cauchy problem ljl.6|l (Section |2J. 

We will use the following notation throughout this paper. 

Notation. If (a £ ) £ g]o,i] an d (& e ) e e]o,i] are two families of numbers, we write 

a £ < b £ 

if there exists C independent of e G]0, 1] such that for any e g]0, 1], a £ < Cb £ . 

2. The Cauchy problem 

Before studying semi-classical limits, we recall some known facts about the initial 
value problem (|l.tjfl . We will always assume that the initial datum / is in the space 
E defined by 

E := {4> G H x (R n ) ; ||0|| s := \\<P\\ L . + \\x^ + ||V0|| L » < +co} . 

This space is natural in the case of Schrodinger equations with harmonic potential, 
since E is the domain of \J — A + |a;| 2 (see for instance [37] ). Local existence re- 
sults for (|1.6(l follow for instance from Strichartz inequalities (one can do without 
these inequalities, see |2Zj)- Global existence results then stem from conservation 
laws (see l|2.3|l below). From Mehler's formula (|1.8fl . Strichartz type estimates are 
available for 

e - 4 i(- £ 2 A+^) = . W{ ^ 

Definition. Let n > 2. A pair (q,r) is admissible if 2 < r < (resp. 2 < r < 
oo if n = 2) and 

- = 5{r) =n - - - 
q \2 r 

Following [jy, we have the following scaled Strichartz inequalities: 

Proposition 2.1. Let I be a finite time interval. 

(1) For any admissible pair (q,r), there exists C r (I) such that 

(2.1) e^\\W(t)cf>\\ Lq{I . Lr) <C r (I)\\ct>\\ L2 . 

(2) For any admissible pairs (qi,r±) and (92,^2)7 there exists C ri)I - 2 (7) such that 



(2.2) e«+^ 



U £ (t- s)F{s)ds 
in{s<t} 



< CVi.raCO W F \\l"2{I-V J 2) ■ 

L1l(I;L r l) 



The above constants are independent ofe. 

The main result of this section follows from ^] . Denote 

Y(I) = {<(> G C(I, E) ; 0, \x\4>, V x cf> G L« oc (7, LQ, V(«, r) admissible} . 
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Proposition 2.2. Fix e e]0, 1] and let f E £. Then (|1.6|l has a unique solution 
u e € y(M). Moreover, the following quantities are independent of time: 

Mass: \\u £ (t)\\ L 2 , 

(2 - 3) i5n e r 32/ : i||eV^(t)||i 2 + i||^W||| 2+£ / (M^ 7 * \u s \ 2 ) \u £ (t,x)\ 2 dx 

It was noticed in |H] that this result can be retrieved very simply thanks to the 
following lemma, which we will use to prove asymptotics. 

Lemma 2.3 ( 6 ). Define the operators 
x 

(2.4) J £ (t) — — sin t — i cos £\7 X ; H £ (t) = xcost + iesmtS/ x . 

e 

J £ and H £ satisfy the following properties. 

• They are Heisenberg observables: 

(2.5) J £ (t) = -iU £ {t)V X U £ {-t) ; H £ (t) =U £ {t)xU £ {-t) . 

• The commutation relation: 



(2.6) 



x? 



H £ (t),ied t + -A-^- 



= 0. 



• Denote M £ (t) = e ^f7 tant ; and Q £ (t) = e i: k cott , then 

(2.7) J £ (t) = -icostM £ {t)V x M £ (-t) ; H £ (t) = ie sint Q £ (t)V X Q £ (-t) . 

• The modified Sobolev inequalities. Let 2 < r < (2 < r < oo if n = 2); the 
exists C r independent of e such that, for any (f> E £, 

ii^iu^aicosr^ii^ii^ii^W)!!^, 



\ L r <a|esini|- 4 W||^||^ S(r) ||^(t)0|| 4 - ( - r) 



(2.. 

ii. . ,, ,,,,, £ , 

• Action on nonlinear Hartree term: for (j) — (f>{t,x), 

(2.9) J £ (t) ((\x\-i * |^| 2 ) 0) = (\x\-> * |0| 2 ) J e (i)0 + 2 Re * (^J £ (i)0)) . 

TTie same /io/ds for H £ (t). 

Remark. Property l|2.6l) follows from (|2.5I) . which is the way J E and ff e appear 
in the linear theory (see e.g. [301 P- 108]). Property i|2.8fl is a consequence of 
Gagliardo-Nirenberg inequalities and l|2.7[) . Finally, l|2.9(l stems from (|2.7|l . 

3. "Very weak nonlinearity" case 

In this section, we study the semi-classical limit of u £ when 7=1 and a > 1 
which is equivalent to "very small" data in our context (cf. (11.21) 1. This case 
includes the 3D Schrodinger-Poisson equation with "very small data" . We prove 
that the Hartree term plays no role at leading order. 

Proposition 3.1. Let f S T,, n > 2, and assume a > 7 = 1. Then for any T > 0. 

IK - wLelli»°([0,T];i=) = ( £ "~ lln 7) aS £ ~* °' 

and for any 8 > (a — 1 — 6 > 0), 

\\A £ (t) (u £ - ulJ\\ Lco{[0>T] . L2) = (e"- 1 - 5 ) as e -> 0, 

where A £ is either of the operators J £ or H £ , and uf lee is the solution of (|1.7f) . 

Remark. Using modified Sobolev inequalities l|2.8|) , we can deduce L p estimates for 
u £ — M| ree for 2 < p < 2nj (n — 2) (2 < p < 00 if n = 2) from the above result. 
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Remark. We could probably get the logarithmic estimate for the second part of 
the statement as well, using Strichartz estimates. The proof given below is not 
technically involved, and suffices for our purpose: we do not seek sharp results. 

Proof. Denote w £ — u £ — u| rcc . ft solves the initial value problem 

ied t w e + ^e 2 Aw £ = ^-w £ + e a (\x\^ * \u £ \ 2 ) u £ ; wf t=0 = . 

Standard energy estimates for Schrodinger equations yield 
(3.1) ed t \\w s (t)\\ L2 <e a \\(\x\- 1 *\u s \ 2 ) 
From Holder's inequality, we have 

1 1 _ 1 

\ L 2 - m*i - i" i \\ L r II" ' ^ r + k~ 2 
From the Hardy-Littlewood-Sobolev inequality, 



u l| L2 



(3.2) \U\x\- 1 * \u £ \ 2 ) u £ \\ L2 < \\lx\- 1 * K| 2 |L. \\u°\\ Lk , for - 



(3.3) \Wx\- 1 * |u e (t)r rr < \\u e {t)\\i P , for Kr,~ < oo and 1 + - = - + -. 
n iiiv - 2 r p n 

Therefore, 1)3.1)1 yields 

(3.4) ed t \\w £ (t)\\ L ,<e a \\u £ (t)\\ 2 LP \\u £ (t)\\ Lk , 

where p and k satisfy the properties stated in (|3.2|l and 1)3.3)1 . For k = 2, r = oo 
and p = 2n/(n — 1), the algebraic identities stated in l|3.2|) and (|3.3|) are satisfied. 
Now since the conditions 1 < r < oo and 1 < p/2 < oo are open, a continuity 
argument shows that we can find p and k satisfying all the properties stated in 
()3.2|) and 1)3. 3JI . Notice that they imply the relation 25(p) + S(k) = 1, hence S(p), 
S(k) < 1; this allows us to use weighted Gagliardo-Nirenberg inequalities. 

We have w £ t=0 = 0, and from Proposition 12. 21 w £ £ C(R + ; £). Therefore, there 
exists t £ > such that 

(3.5) \\J £ (t)w £ \\ L 2 < 1, 

for < t < t £ . The argument of the proof then follows [2E] (see also [H]). Recall 
that from Lemma IP || J £ (t)u e ilee \\ L 2 = ||V/|| L 2. 

Because of l]2.6p . J e u| rcc solves the linear Schrodinger equation with harmonic 
potential, and ||</ e (£)uf roc |U 2 = ||V/||l2. So long as (|3.5|l holds, we have, from 

ll^)ll-<^y ; ll«m.<|^, 

for some Co independent of e and i. Then 1)3. 4|l yields 



e9 t ||«; £ W|| i2 



< 



\coBt\ M <P)+*W |cost| 
Integration in time on [0, t] yields, so long as 1)3.5)1 holds, 

" 0lT 



\W ||L~([0,t];L 2 ) 



< 



a-1 



COST 



For t < it/2, we get, so long as 1)3.5)1 holds, 

ll^||L~([0,t];L 2 ) ^ Q_1 

From 1)2.6)1 . J £ (t)w £ solves 



M§-« 



ied t J £ w £ + ]-e 2 AJ £ w £ = M- J E w £ + £ Q J £ ((M" 1 * K| 2 ) m £ ) ; J E wf t=0 = . 
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Using i|2.9jl. energy estimate for J £ w £ yields 
ed t \\ J e (t)w £ \\ L 2 < e a (\K\x\- 1 * K| 2 ) J £ (t)u £ \\ L2 + \\lx\- 1 * (WJ £ u £ ) ■ u £ \\ L2 ) 

< e a (lllxp 1 * K| 2 || Loo \\J e {t)u £ \\ L , + \\lx\- 1 * (Wj £ u £ ) • u e || L2 ) 
For the first term of the right hand side, use the easy estimate 

||N- 1 */||^II/IIl(»-)' + II/II 1 c»+)' 

where n~ (res. n + ) stands for n — r\ (resp. n + rj) for any small rj > 0. We have 

9n 

IIN- 1 *Ki 2 l| i0O <ik e wn 2 Lre - + iKwiii, + , with k = — x 



It is at this stage that we lose the logarithmic rate (we cannot use Hardy-Littlewood- 
Sobolev inequality when an exponent is infinite): using Strichartz estimates (see 
Section EJ, we believe that we could recover that rate, with a more technically 
involved proof. 

For the second term, we proceed as in the beginning of the proof. From Holder's 
inequality, 

(3.6) *lFJ £ u £ ) ■ u £ \\ r2 < II Id" 1 * (TiFJ E u s )\\ TT \\u £ \\ L , , with - + - = -• 
" / in* ii ii X/ r a 2 

From the Hardy-Littlewood-Sobolev inequality, this is estimated, up to a constant, 
by 

(3.7) \\WJ £ u £ \\ rp \\u £ \\ L <, , with 1 + - = - + - forl<r,p<oo. 

r p n 

Using Holder's inequality again yields an estimate by 

1 _ 1 1 
p ~ 2 + k ' 

Take r = n, a = 2n/(n — 2), k — 2 and p = 1: the algebraic identities from f>|) . 
(|3.7|) and l|3.8l) arc satisfied, but not the bound p > 1. Decreasing slightly a increases 
p (take a large but finite when n — 2), so we can find indices satisfying (|3.6() . (|3.7() 
and (|3.8I) by a continuity argument. Note that they satisfy 8{k) + 5(a) = 1, and 
each term is positive. 

Gathering all these estimates together yields the energy estimate 

Ed t \\J £ (t)w £ \\ L , < e a (\\u £ (t)\\l^ + \\u £ (t)\\l K+ + \\u*\\ Lh \\u £ \\ L „) \\J £ (t)u £ 

So long as l|3.5|l holds, we deduce from i|2.8[l . 

1 1 1 



(3.8) ||u e || ifc || J £ u e \\ L2 \\u £ \\ L * , with 



L 2 



ed t \\J £ (t)w £ \\ L 2<e a 



| cost| 2<5 ( K "> I cost| 2 < 5 (« + ) |cost|' 5 ( fc )+ 5 ( ff ) 



• -" ( 4- — L_ \ < 



\cost\ 2S ^ + ) I cosily ~ |cosi| 1+ ' 
Integrate this, so long as (|3.5(l holds: 

ll^llL-([o J t] ii =")<e a - 1 (|-*)° 
Fix 5, A > 0. So long as (j3.5jl holds, we infer, for t < n/2 - As, 

\\J £ w £ \\ L oo m] . L 2 ) <e a - 1 (AE)- S . 

Therefore, there exists ea > such that, for < s < ea, l|3.5|l holds up to time 
7r/2 — Ae, with the estimates 

(3.9) ||w £ || i oo ( [ 0Tr/2 _ A£ ]. i 2 ) < £ Q_1 In - ; || J £ W £ \\ L °a([ ,Tr/2-Ae];L 2 ) S E a ~ 1 ^ 5 ■ 
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An estimate similar to that of J £ w £ then follows for H £ w e , since from 12.3fl . 
\\H £ (t)u £ \\ L * < 

Denote I £ K = [n/2 — Ae, ir/2 + Ae]. Mimicking the above computations, we have 



\\w £ \\ L ^(i' A -L 2 ) 



< 



77 A 

2- Ke 



L 2 



\u%T)\\U\u £ (T)\\ Lk dT, 



where p and k satisfy (|3.2[1 and i|3.3[l . Recall that they satisfy 2S(p) + S(k) = 1. 
Using the conservations of mass and energy l|2.3[l , along with Gagliardo-Nirenberg 
inequalities, we have, for any t, 

\\u e (t)\\ LP <s- s ^ ; \\u%t)\\ Lk <e-^ k K 



We deduce 



\ w \\L^{Il;L 2 ) 



< 



ur 



( - - Ae) + £ ^ £ -^(P)S(k) ,je | < £ «-i ln I + Ae o-i 
V 2 ' l 2 e 



The same method yields, since l|2.3l) shows that ||-fP(t)M e ||z,2 < ||/||s: 

< , for any 5 > . 

To treat the case of J £ w £ , introduce 

z&)= sup ||J £ (rK|| L2 . 

■|-Ae<r<t 

Proceeding as above, we have 



< 



(3.10) 



L 2 



-Ae 



\J £ (t) (l^r 1 * \U £ \ 2 U £ ^ 



\L 2 



(It 



< e a - 1+ +e a - 1 



e- 1 (z £ (r) + \\J £ ( T )ul cc \\ L2 )dT. 



-Ae 



We can then apply the Gronwall lemma (recall that || J e (r)u| roc ||i2 = ||V/||i2): 



< 



Gathering these informations we get, for any d > 0: 



\\J £ w £ \\ LX{Il , L 2 ) + \\H £ w £ \\ Lx ^ L2) <e c 
For t G [n/2 +e,tt], we can use the same proof as for t € [0, ir/2 — e], to obtain: 



\W ||Loo([ 0j7r ] ;L 2) ^ e 



<e" _1 ln-. 



m E ^ S |U°°([0,7r];i 2 ) + \\H e W £ \\ L ^^^]:L 2 ) 



Repeating the same argument a finite number of times covers any given time interval 
[0,T] and completes the proof of Proposition l3.il □ 

4. Nonlinear propagation and linear focus 
In this paragraph, we assume a = 1 and 7 < 1. We define 

dr 



(4.1) 



g(t,x) = -(\xr*\f\ 2 )(x) 



costT 



This function is well defined for any t, since 7 < 1. We will see later on how this 
function appears. 
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Proposition 4.1. Let n > 2, / e S, and assume 7 < a = 1. Let A e be either of 
the operators Id, J e or H € . 

• For < t < tt/2, the following asymptotics holds: 



sup 

0<T<t 



A' 



(r) ( u £ (t, x) - 1 ——f ( —) e"*^ tanT +^( T ^) N j 

w V (cost)"/ 2 ^ Vcosr/ / 



•0. 



For tt/2 < t < tt, 



sup 

t<T<7l 



A' 



w V (cost)"/ 2 ^ Vcosr/ 



L 2 „ 



6^0 



For t = ir/2, 



B e u 



(D-X^ W )G) 



L 2 



£^0 



0. 



where B e is either of the operators Id, - or eV x , and the Fourier transform is 
defined by 

(4.2) m0 =m = tA^H I e~ lX< H*)dx ■ 

\^) 1 Jmp 

Remark. We can also prove estimates for arbitrarily large time intervals, with the 
same proof as below. 

Remark. The difference between the asymptotics before and after the focus is mea- 
sured only by the Maslov index, through the phase shift e - m7r / 2 : n o nonlinear 
phenomenon occurs at leading order near the focus. On the other hand, nonlinear 
effects arc relevant outside the focus, as shown by the presence of g. 

4.1. Oscillatory integrals. The main tool for proving Proposition l4. H is the same 
as in linear cases ([3], see also [21 EH for applications in nonlinear settings): we 
represent the solution u e as an oscillatory integral. Recall that u £ € C(R; S) and 
that e _I A( _e A+x ) = U e {t) is a unitary group on L 2 . Define a £ by 

(4.3) a e {t,x) = U £ (-t)u e (t,x) . 

We first seek a limit as e — > for a £ before the focus. This is suggested by a formal 
computation as in [3], and the following lemma: 

Lemma 4.2. For t £ [0, 7r/2[U]7r/2, tt], define V s by 



(4.4) V s (i)0(a;) 



(cost) 



cost 



tW 2<?!> (cos t) 



™/ 2 9 V cos < / 



-VWK 1% 



i/0 < * < tt/2, 



if ir/2 <t <n . 



For any (/) € H 1 (R n ), any 6 e]0, 1/2], and any i e [0, 7r/2[U]vr/2, tt], 

||W £ (t)0-V E (i)0|| i2 <2|etant| e || ( /)|| /fl . 
Proof. Notice that from Mehler's formula 1(1.8(1 . we can write, for < t < tt, 



W(t) = MPZFMI where M e t {x) = 



D^(x) 



1 



{iesmt) n / 2 Vsint/ ' 



and the Fourier transform is defined by 14.2(1 . We infer 



\\w (t)cf> -r m 



L- 



1 



From Parseval formula 
1 



(2z7rtani)™/ 2 



(2Mrtant)V 2 

■ |x-y| 2 



1 



(2tt)™/ 2 



e H*m*f[ v )dy -f(x) 

- ietant ^r +lx <Tf{i)d^. 
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therefore 



(2tt)™/ 2 



-is tan t %- 



*■ - 1 J7(0 



L 2 



from Plancherel formula. The lemma then follows from the estimate \e ls — 1| < 2\s\ e , 
forO<(9<l/2. □ 

From Duhamcl's principle, we have 

u £ {t) = U £ {t)f - i f U £ {t - s) (OP * \u £ \ 2 ) u £ ) (s)ds . 
Jo 

Using i|4.3fl. we deduce 

(4.5) d t a £ (t) = -iW(-i) * K| 2 ) if) (t) . 

Now the formal computation begins. Assume a £ — + a as e — * 0, in some suitable 
sense. Then u £ (t) ~ U £ (t)a(t), and from Lemma f4. 21 



u £ (t,x) ~ — ft, - ) 

V (cost)"/ 2 V COSt/ 



e _ lf _ tant for < t<7r / 2 



Plugging this into (|4.5(l and using Lemma ^21 again (with U £ (—t) instead of U £ (t), 
the result still holds) yields 



d t a(t,x) 



costp 



(\xr*\a\ 2 )a(t,x). 



Recall that a\ t =o = uf t _ — /, and notice that from the above ordinary differential 
equation, 9t|a| 2 = 0: we have a(t,x) — f(x)e , ' 9 ^ t ' x \ where 



d t g(t, x) 



■1 



COStT 



Or* |/l 2 ) (*) ; sit=o = o. 



Integrating this equation yields the definition of g(t,x) given in l|4.1|l . 

Proposition 14.11 stems from the more precise following proposition, Lemma 14.21 
and a density argument. In view of a rigorous justification, denote 

(4.6) b £ (t, x) = a £ {t, x)e~ i9 ^ = e^'^W : (-t)u e : (t, x) . 

Proposition 4.3. Let f E £ n H 2 (R n ). Fix 5 > 0. There exists C s such that 

sup \\b £ (t)-f\\x < [ \\d t b £ {t)Udt < Cse 1 -'- 5 . 

0<t<7T JO 

The first inequality is trivial. We prove the second one in three steps: 

(i) On [0, 7r/2 — Ae] for any A > 0, with a constant depending on <5 and A. 

(ii) On [tt/2 — Ae, tt/2 + Ae], with a constant depending on S and A. 

(iii) On [tt/2 + Ae, tt], with a constant depending on 8 and A. 

As in Section [31 the parameter A > is arbitrary, while it has to be large in the 
case a = 7 > 1 (see Section [S] and [S]): this situation is typical from a case where 
the focus is "linear". 
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4.2. Asymptotics before the focus. Fix A, 6 > 0. We prove that there exists 
Ca a such that 



r 



-Ae 



(4.7) 
Denote 

Ve{t) = 

From (|4.5(l and the definition (|4.6(l 



d t b £ (t)\\vdt<C A ,s£ 



l-7-<5 



\d t ¥{T)\\ m dT. 



\\dtb e {t)\\ L , = 
(4.8) 



W{-t)((\x\-~<*\u^)uZ){t)- 



(h-^i/ivw 



(M 



j cosi|' 



, cosi| 7 
W(t)((\xr*\f\ 2 )a s ) (t) 



L- 



V 1 



Lemma 14.21 suggests that we can replace U e with V £ in the last expression, up to 
a controllable error. Before going further into details, we prove two lemmas which 
will be of constant use in the proof of Proposition 14. HI 

Lemma 4.4. Assume 7 < 1, and let < 6 < 2(1 — 7). There exist p and q with 

n 



2S(2p')= 7 +°-{<l) ,p<- 
2 7 



W) = i±i + |(<i), 9 < 7+1 , 



and such that there exists C such that for any <f> £ C£°(R n ), 

||M- 7 *0|| LOO <c(U\\ L i + u\\ Lpl ) , 
||v(M- 7 *^)|U <c(\\4>\\li + \\4>\\l«>) ■ 

Proof. We have 2S(2p') — 7 when p — n/7, and S(2q') = when q — n/(j + 1). 
Therefore p < n/7 and g < 71/(7 + 1) if 26{2p') =7 + 5/2 and S(2q') = 2+i + |. 
Let x 6 C~(M+, [0, 1]) with x = 1 on [0, 1]. We have 

111*1-7 < ||(xkr 7 ) *0|| L « + ||((i-x)kr 7 ) 

< llxN -7 !!^ 11011^ + 11(1 -x)N- 7 IU-H0IUi 



where we have used x i— > \x\ 
similar, since V|x| _7 = &{\x\ 



£ Lf C (K W ) because p < n/7. The other estimate is 



□ 



Lemma 4.5. Le£ 7 < 1 and f e Sni? 2 (M n ). iJecaZi that g is defined by (|4.1|l . We 



M" 7 * |/| 2 G VF^ 00 ; ff S L£ C (R; W 2 ^) ; /e", (M' 7 * |/| 2 ) /e'» G L£ C (R; # 2 ) . 

Proof. From Lemma 14 . 41 and Sobolev embeddings, 

2 

\L 2 p' 



\xr*\f\ 2 \\ L ^<\\f\\h + \\f\\ 2 L2P ,<\\f\\ 2 m , 



|V|s 
721 



L-Sll/lli- 



< 11/11^ 



1 1 ■ * 1/1 2 _ 

iv^krH/I'Loo^llvw-^^vi/i 2 )^ 



Since t 1— ► | cos £| 7 e Lio C ( 
properties follow easily. 



<ll/IU»l|v/|U a + ||/|| £ v||v/|| za ,/ <||/||| 2 . 

we infer that g € ££ C (R; VF 2 * 00 ). The last two 



□ 
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We can now replace U £ with V e in Q4.8[l . up to the following error. From Lem- 
mas Ol El and E31 



\{W(t) - V £ (t)) ((\x\-i * |/| 2 ) a e ) (t)|| i2 < |etant 

< |etant 

< jetani 

< |etant 

< jetani 

< |etant 

< |e tani 



N" 7 *l/| 2 ) a e (t)\ 



ff 1 



l/l 2 



K(*)ll 



(||a £ (<)|| i2 + ||V^(i)|| L2 ) 
(||/||L» + ||V a (6 e ^)|| i2 ) 
(l + HV^WHiO 
(1 + ||V, (b%t)-f)\\ L2 ) 



1+ ll^& £ (r)|| ffl dr , 



for < 8 < 1/2 to be fixed later. Plugging this estimate into 14.8JI yields 

1 



II WO II L 2 



< 



(4.9) 



x\-~< *\u e \ 2 )u e 



(t) 



j cost| 



;V £ (*)((M- 7 *I/I 2 ) « £ ) (*) 



L 2 



|etani| 6 
|cost|T 



We check that 



(4.10) 



| cost|' 



-Xii) {{\x\-i * l/l 2 ) 0) = (M~ 7 * |v e W/| 2 ) r (t)* 



Since we expect V £ (t)a £ (t) to be close to U £ (t)a £ (t) = u £ {t) as e — ► 0, we estimate 
the difference 



(\x\-y*\r(t)f\ 2 )(nt)a £ (t)-w(t)a s (t)) 
<|||zp*|v £ (t)/| 211 



L 2 



ff 1 



L ~\\(V e (t)-U £ (t))(b £ e^)\\ L2 

< (\\^(t)f\\h + l|V e (t)/||iv) (^an^ |&^)e^> 

< (l + | cosi]- 25 ^ ) ktani| 9 (\\b £ (t) - /|| H1 + H/lUx) , 

using the modified Sobolev inequality l|2.8[) . Since 2S(2p') = ~ > 7, we infer from 
jOjl that 

(4.11) 



iw*)iu*< , |g ^L a+fe(*)) 



|cos^| 2 W) 

+ ||(N-' lf *(|« e (t)| 2 -|v-(t)/| 2 ))« e (t)|| La ■ 

From Lemma 14.41 the last term is estimated, up to a constant, by 

||K(i)| 2 - |V £ (i)/| 2 || L1 + \\\u £ (t)\ 2 - |V^)/| 2 || LP , < 
(4.12) £ ||« £ W - V £ W (/e iffW ) || i2 (||« e (*)|| £ , + |v £ W (/e*«) || l2 ) 

+ ||« e (t) - V £ (t) (/e l9(t) )|| i2p , (ll^WIUv + ||v £ W (/^ (t) )|| i2p ,) 
For the first term of the right hand side, we have, since IA £ is unitary on L 2 , 
\\U £ {t) (b £ e^) - V £ (t) (/e*») || i2 < \\b £ (t) - f\\# + \\(U £ (t) - M £ (t)) (je*®) 



L 2 



< 



y £ (t) + |etant| £ 



fe 



ig(t) 



i? 1 
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In addition, notice that ||M e (t)|| L 2 = |jV e (i)/|| L 2 = ||/||j,2. The second term is 
estimated thanks to the modified Gagliardo-Nirenberg inequality (|2.8(l , 



(t) - V £ {t) (/e^W) I < | cost\-^ 2p ^ ju £ (t) - M £ {t) (/e l9(t) ) 

J £ (t) (u e (t)-V s (t) (/e i9 «)) 



l-5(2p') 



L 2 
S(2p') 

L 2 



The first L 2 -norm was estimated just above. For the second one, notice that 



J £ {t)U £ (t) = -iU £ {t)V x ; J £ (t)V £ (t) = -i V £ (t)V x , 



therefore: 



:(t)(u £ (t)-V £ (t)(fe^)) 



< 



W 



< 



(t)V (V(i)e i9(t) ) - V £ (t)V (/e l9(t) 
(U £ (t) - V £ (t))V (fe t9{t] 



L 2 



v(V(i)e i9(t) -/e l9(t) ) 
< j/ £ (i) + |etani| e , 
where we have used Lemmas 14.21 and 14.51 We infer that 

\u £ (t)~V £ (t) (/e l9(t) )| i2p , < |cost|- 5 ( 2p,) (y £ (t) + \ettmt\ e ) . 
We have explicitly 

V £ W (/e*<'>) , =|cost|-W|| / || < Icos^r 5 ^'). 

V / L 2 p 

Proceeding as above, we have 

\\um)\\ L2p , < |cosr 5(V) ii« £ ii^ (V) iu £ W" £ ii^' ) , 

with || J e (t)M e ||i2 < ||6 E (t) — /Wh 1 + ll/llff 1 - These estimates will eventually lead to 
an inequality of the form y' E (t) < a(t)y e (t) + b(t)y e (t) K + c(i), for some « > 1. To 
avoid that situation, we proceed as in Section [3| there exists t £ > such that 

(4-13) \\b £ (t)\\ H i < 211/11*1, 

for t S [0,t £ ]. So long as (pTTgj) holds, we have from the above estimates 

(4.14) \\d t b £ (t)\\ L 2 < | cost|- 25(V) (y £ (t) + | e tant| e ) . 

To prove that Ij4.13|l holds up to time tt/2 — Ae for < e < £\ along with the error 
estimate l[4.7jl . we estimate the L 2 -norm of W x dtb £ . From (|4.5(1 and l|4.6|l . 

V x dtb £ (t) = -iV x g(t)d t b £ (t) 

1 



ie~^V x W(-t) ((\x\-i * \u £ \ 2 ) u £ ) (t) 



I COS^|T 



{\x\-^\f\ 2 )a £ (t) 



The first term is controlled thanks to Lemma 14.51 and (|4.14|) . For the other term, 
we notice that since U £ is unitary on L 2 , from (|2.5|) its L 2 -norm is equal to: 



J £ (t)((\x\-^\u £ \ 2 ) u £ ) (t) + 



cosip 



U £ (t)V x {{\xr*\f\ 2 ) a £ (t)) 



We proceed as before: we first replace IA £ with V e in the last term, up to an error 
of | cost| 1 times: 



(U £ (t)-V £ (t))V x ((\xr*\f\ 2 ) 



< 



L 2 



< 



\{U £ {t) - V £ (t)) X7 X ((|x|" 7 * |/| 2 ) (b £ - f)e ig ) | 



+ \\(W(t)-V £ (t))V x ((\xr*\f\ 2 )fe i *) 
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For the first term, we do not use Lemma f4. 21 but roughly the fact that U e and V e 
are unitary on L 2 . It is not larger than 

2 ||V, ((|:rr * l/l 2 ) (b £ ~ f)e ig ) \\ L2 < W(t) - f\\ m , 

from Lemma f4. 51 The second term is controlled thanks to Lemmas 14 . 21 and 14.51 

|| (W(t) - V%t)) V, ((\xr * l/l 2 ) fe lg ) \\ L2 < \e tani| e . 

We now have, so long as (|4.13|) holds, 



iiwoii 



< 



(4.15) 



it|- 2W) (ys(t) + \etant\ e 



j £ (<)((N^*KI 2 )« e ) + 

Using the identity J e (t)V £ (t) 

r(t)({(\ x r*\u*\ 2 ) u s ) 

(4.16) 



r(t)V x {(\x\-i*\f\ 2 )a*(t)) 



,COSt|T 

iV £ (f)Va; and (|4.1U|) . we have to estimate 
1 



,COSt|T 

r(t) ((\xr * ki 2 ) u « - (\x\-i * \r(t)f\ 2 ) r(ty) 



L 2 



L 1 



\L 2 



< || {\xr * \u^\ 2 ) jw - (m- 7 * iv £ w/i 2 ) j £ (t)v e (tK|| i2 

+ | cosi| ||V a (\x\-i * K| 2 ) - V x (\x\^ * \r{t)f\ 2 ) r(t)a e \\ L2 . 
We replace V e with hi 6 in the first term of the right hand side, up to the error 
|| (|*r * |V E (*)/| 2 ) J 5 {t) (V £ (i) - Wit)) a^\\ L2 < 

<(PHt)f\\h + \\r(t)f\\ 2 L2p ,)\\(r(t)~u%t))v x a^\\ L2 

<\ cost\- 2S ^ || (r(t) - W{t)) V x {{¥ - fie*) \\ L2 
+ | cost\- 2S ^ \\(r(t) - W{t)) V x ((/ e *») || i2 

<icosr 2W) (n &e (*)-/ii^ + i £taii *i e ) . 

from the above computation. Therefore, the first term of the right hand side of 
(|4.1(i|) is estimated by 

| cost|- 25(2p,) (lfe(«) + \etaat\ B ) + \\(\x\^ * (K| 2 - \V E (t)f\ 2 )) J E (t)u s \\ L2 . 
So long as Q4.13|l holds, || J e {t)u e \\i J 2 < 1, and the last term is estimated by 



(h 



|v E (t)/| : 



which already appeared above and was estimated in H4.12JI . We are left with the 
second term of the right hand side of l|4.16[) . Using Lemma T4 . 41 with q instead of p 
now, 



|v* (i^r * \r(t)f\ 2 ) (v e (*) -w{t)) a z\\ L2 
< 



< 



(\\^(t)f\\ 2 L2 + \^(t)f\\ 2 L2ql )\\(r(t)^W(t))a^\\ L2 
<\cost\- 2 ^ 2q '^ (y £ {t) + \et£mt\ e ) . 
The final term to estimate is 

||v x (M~ 7 * (Kl 2 - |v £ (t)/| 2 )) u s \\ L2 < ||K| 2 - \r(t)f\ 2 " 



IL 1 



,£|2 



|v £ W/| 2 



\Li' ■ 



The right hand side was already estimated in (|4.12(l with p instead of q. We finally 
have, so long as (|4.13jl holds, 

y 6 (t) + \et&nt\ ) . 
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Now recall that given 6 > 0, 8(2p') and 6(2q') are explicit, hence 
y\t) < |cost|-T-4 (y e (t) + |etani| e ) . 

It is now time to fix 9. In view of l|4.7Jl . it is natural to take 8 — 1 — 7 — 5. This 
yields, so long as H4.13|l holds, 

(4.17) y' £ {t) < |cosir~ f {y E (t) + \et^t\^- s ) < j cost|^-| y E (t)+ , ■ 

I cos 1 1 ^ 2 

The maps £ 1— > |cosi|~ 7- ^ and i 1— > |cost|~ 1+ 2 are locally integrable (we can 
assume 7 + 6/2 < 1 — 5/2, otherwise Ij4.7|l is of no interest). From the Gronwall 
lemma, so long as (|4.13|) holds, we infer 

(4.18) ysW^e 1 -*- 5 . 

Therefore, there exists £a > such that for < e < £a, <|4.13|1 holds up to time 
7r/2 — Ae, with Ij4.18|l . The estimate for xdtb £ then is easy, we leave out this part; 
this proves (|4.7jl . 

Remark. One might believe that we could deduce Proposition ^. 31 in one shot from 
(|4.17l) . and wonder why we split the proof into three steps. The reason is that we 
cannot apply Lemma l4~2*l f which was used to get Ij4.17|l ) near t = ir/2. On the other 
hand, we will see below that computations near t = tt/2 are far simpler. 

4.3. Near the focus and beyond. Keep A, S > fixed. We prove that there 
exists Ca,s such that 

(4.19) [ 2+ " \\d t b £ (t)\\zdt < Cajs 1 -'- 5 . 



11 Wi)ik2 < || (\ x r * ki 2 ) u £ (t)\\ L2 + — — 11 (\ x r * 1/1 2 ) « £ w|| L 2 



A rough estimate in l|4.8|l yields 

1 

IcosilT 

(4.20) 1 1 

< (\\u s (t)\\h + \\u%t)f L2p ,) \\u%t)\\ L * + \\u £ (t)\\ L2 . 

The conservation of mass yields ||it £ (i)||£2 = ||/||z,2. The conservations of mass and 
energy Ij2.3[l yield, along with Gagliardo-Nirenberg inequalities, 

\\u e (t)\\ L2pl <e- s ^'K 

Using this estimate (which is sharp near the focus, and only near the focus) and 
integrating H4.20J) . we get 

T+Ae ri+ Ae rlf 

\\d t b e (t)\\ L2 dt < ke 1 - 2 ^ + 1 ' 



Ae " J%-Ae \COStp 

The term ||iE<9 t 6 £ (t)||£2 is estimated the same way, since the conservation of energy 
yields an a priori bound for H e u e . For || V x dtb s (t)\\ li , we proceed as in Sectional 
(|3.1U|) to get an estimate from Gronwall lemma; the details are left to the reader. 

Finally, one can prove that there exists Ca,8 such that 

Wdtb'm^dtKCXse 1 ^- 6 



l+Ae 

by mimicking the computations performed in Section l4.2l and the proof of Propo- 
sition is complete. 
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5. Linear propagation and nonlinear focus 

We now consider the case where a = 7 > 1 in 1)1.6(1 . Our results are similar to 
those of pjj. Before stating the main result, we recall some points of the scattering 
theory for (|1.9|) . 

Proposition 5.1 ( [I'M I17p . Assume ip- € E and 1 < 7 < min(4, n). If 7 > 4/3 
or i/UV'-lls *s sufficiently small, then: 

• There exists a unique ip G C(Rt,E) solution of <|1.9|) . suc/i £/iai 

lim IIV- -U(-t)-0(t)||s = 0, where U(t) = e^ A . 

t—> — oc 

• There exists a unique ip+ G S swc/i 

u(-t)v(t)|| E = o. 

TTie scattering operator is S : ?/>_ i— > '0+ ■ 
Our main result in this section is: 

Proposition 5.2. Suppose n > 2. Let f € E, 1 < 7 = a < min(4, n), and fc G N. 

Assume either 7 > 4/3 or ||/||s *s sufficiently small. Then the asymptotics of u € 
for tt/2 + (k — l)7r < a < b < tt/2 + kit is given by 



lim ||V>+ 
/ — '+00 



sup 

a<t<b 



A £ {t) u e {t,x) 



-inki 



it\ n / 2 



(ToS k or 



' \ cos t J 



L'i 



0. 



where A e is either of the operators Id, J 6 or H e , and S denotes the k-th iterate 
of S (which is well defined under our assumptions on f). At the focuses: 



+ kn) - 



1/2 



(ToS k )f 



L 2 



0. 



where B e is either of the operators Id, - or eV 2 



With Lemma |4.2I in mind, this shows that nonlinear effects are negligible away 
from focuses, while they have an influence at leading order near the focuses: each 
caustic crossing is described in average by the nonlinear scattering operator S (the 
phase shift e ~ mfe 5 is the Maslov index, present in the linear case [Hj)- 

The proof of Proposition 15.21 is very similar to the one in , which relies on 
(scaled) Strichartz estimates. We will refrain from repeating everything in detail 
and limit ourselves to prove the main technical proposition and present an outline 
for the rest of the proof. One main difference to the problem in |H| is the action of 
the operators J e (t), H e (t) on the Hartree nonlinearity as described by l|2.9|l . 

We start by reformulating Equation l|1.6|l by the Duhamel formula 



u £ {t) =U e (t- t )u e 



(5.1) 



ie 



7-1 



W(t - s)F E (u £ )(s)ds 



U e (t- s)h £ (s)ds. 



This equation generalizes Eq. (|1.6f) to the case of an additional source term and a 
general nonlinear term F e . The main technical result which is used throughout the 
proof of Proposition 15. 21 is: 
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Proposition 5.3. Let t\ > to, with \ti — to\ < it. Let q,r,s,k £ [l,oo] be such 
that: 

..1127 , 2n 

(a) — = — I 1 1 and s < 



(5.2) 



r r s n n — 7 

». 1 12 

(c) (q,r) is an admissible pair , 



(d) < i < <*(*) < 1 . 

Assume that there exists a constant C independent of t and e such that for to < 
t<t x> 

(5.3) ||F £ (</)(t)||^ < C — \\u s (t)Ui > 

(| cosi| + e) w 

anci define 



A £ (t 0l h) 



/t (I cosi| + e) fc,5 ^- ) y 
Then there exists C* independent of e, to and t\ such that for any admissible pair 
(P, a), 

ll« B IU«(to,tii£') <C*e- l '"\\ul\\ L 2 +C q<p e- 1 -*--r\\h s \\ LP , {toM . L „- ) 

+ C*e 2 ( d ^^A^t ,t 1 )\\u^\ Lq{toMlLr) , 

Mostly the following corollary is applied: 

Corollary 5.4. Suppose the assumptions of Prov. \5.!A are satisfied. Assume more- 
over that Ce 2 ^^" ) A e (to, £1) < 1/2, which holds in either of the two cases, 

• < to < t\ < J — Ae, with A > Ao sufficiently large, 

• to,ti £ — Ae, ^ + Ae], with ~ < rj sufficiently small. 

Then 

(5.5) ll« e |U-Cto,tx,-£») < C||u§|| ia +C , g , p e- 1 -p||/» e || £(> / (t0itl . i<r /j. 

To prove Proposition 15. 31 we first prove the following algebraic lemma: 

Lemma 5.5. Let n > 2, and assume 1 < 7 < min(4, n). Then there exist q, r, s,k £ 
[l,oo] satisfying the conditions (|5.2|) . 

Proof. Note that (a) is equivalent to demanding 7/2 = 5(r) + 6(s) and 7/2 > 5(s). 
Case 7 < 2: Suppose 7/2 = 5(s). Then, by the first half of (a) 5(r) = and 
(q, r) — (oo,2) by (c). With k = 2, (b) and (d) are satisfied. Now choose s 
such that 1/2 < 5(s) < 7/2, but close enough to 7/2 for (|5.2|l still to be valid by 
continuity ( for example 5(s) — \ + i (Qj- — i) ). Then l|5.2|) is satisfied. 
Case 7 > 2: In this case take s such that 5(s) — 1, e.g. s — -^5. Up to a continuity 
argument as in the previous case, 5(s) < 1 and Ij5.2|l is satisfied. □ 

Proof of Provosition \5.'A Application of the (scaled) Strichartz estimates fProp. l2~T|) 
to equation (|5.1|l yields 

\\u £ \\mt M;Ln ^Ce-^Kll^ + a^e-M-? \\h e \\ Lpl{toM . Lal) 

+^- 1 -fll^(u e )ll^ (t01tl; ^ ) . 

Then by the assumptions on F e (u e ), after an application of Holder inequality in 
time, the statement follows. □ 
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Proof of Corollary \5.4\ The additional assumption implies that the last term in 
(|5.4() can be absorbed by the left hand side, and we get 

(5.6) ||« B |U«(to,t i; iO^ Cfe " 1/9 IKIl£ a + Cre " 1 ~'~'ll AS ll^(*o,t i: £'')- 
Another application of Strichartz estimates to equation l|5.1|l . with indices (oo,2) 
on the left and (p, a) respectively (q, r) on the right, yields 

KIU~(t ,ti;£ 2 ) <C\\u e \\ L 2 +Ce- 1 -p\\h E \\ LP , {to ti . Ltr , ) 

+ Ce^-HF^)\\ Lil{tQM .^ y 

As before, 

< Cev\\u E \\ Lq{to j 1 . Lr) , 
and the statement now follows from (|5.6|) . □ 

The proof of Proposition 15 . 21 consists of three parts: the propagation before the 
focus, the matching between the two regimes, and proof that near the focus, the 
harmonic potential is negligible. In all parts the main tool to derive the major 
statements is Prop. 15.31 Since the proof is very similar to the one in we do not 
repeat everything in detail but give a detailed proof only for the first part to show 
how the methods of are applied. 

We now show the proof for the propagation before the focus, that is the ap- 
proximation of u e (t) by itf ree (£) for < t < -| — Ae, in the limit A — > +oo. We 
prove: 



lim sup sup 

e^O 0<t<£-Ae 



A e (t) (u e (t,x)-u! Iee (t,x)) 



— » o, 

L 2 A— >+oo 



with A e (t) being either of the operators Id, J £ or H e . 
Define the remainder w e = u e — uf ree . It solves 

ied t w E + ^e 2 Aw e = V{x)w £ + e 1 (l^r 7 * |u e | 2 ) u £ , 
«jf t=o =0. 

From Duhamcl's principle, this can be written as 

(5.7) w E (t) = W{ty - is 1 - 1 f U e (t - s) * \u e \ 2 ) u s {s)ds. 

Jo 

Since ztf ree solves the linear equation (|1 . T|) . so does J e (t)uf TCC from (|2.6|l . and 

hl Iee (t)\\ L * = \\f\\» ; II J e (tX b j\ L 2 = ||v/|U 2 . 

From the Sobolev inequality l|2.8(l . 

IKeeWIL* < r^IW ll/II^IIV/ll^ 

for any s G [2, -^jM. Therefore there exists Co such that 

Co 



(5.8) IK-WL. < 



cos 



From Prop. 12.21 for fixed e > 0, u E £ C(R, £), and the same obviously holds for 
u frcc- Therefore, w e £ C(M, £), and there exists f e > such that 

(5-9) IKWllj. < , C °. x , „ 
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for any t 6 [0, t s ]. So long as (|5.9|l holds, we have 

2C 



u 



and we can apply Prop. l531 

T&keh 6 = e 7 (|ir|- 7 * |w e | 2 ) u| eB andF e (w e ) = * |u £ | 2 ) w e and let g, k, r, s <= 

[1, oo] satisfy the assumptions of Prop. 1531 Now by Holder's inequality 

||F £ (0(*)|| ir , < \\\x\^* \u%t)\ 2 \\ LfS \\w%t)\\ L r 

with f3 such that ^7 = 7 + ^ ■ By the Hardy-Littlewood-Sobolev inequality and the 
above estimate, 

lir-KJWII^ < ||«'(i)||i.|K(i)IU- 

Note that the second statement of 1)5.21) fa) ensures that s, (3 € (1, 00) so the Hardy- 
Littlewood-Sobolev inequality is applicable here. Assume l|5.9|l holds for < t < T £ . 
IfO<i<T e <-| — Ae, then e < cost, and the above estimate shows that F e 
satisfies assumption 15.311 . 

From Corollarv l5.4l if A is sufficiently large, we get for < t < T £ < J — Ae: 

IKIU-(o,T;jy») < C^e 1 - 1 '^ \\{\x\-~< * \u £ \ 2 )uf ICC \\ LP , {0 T . La , } 

for any admissible (p,cr). Now take (p, a) = (q,r) and proceed as above in space, 
and apply Holder inequality in time: 

IKl^r 7 * Kl 2 ) U froc|| L ,' (0: T;L'-') - C 7,nlK II i*(0,T;£») IKreelU^O.T;.^) • 

The first term of the right-hand side is estimated through l|5.8|l and l|5.9|l : 

II e 1 1 2 ^ C 



L*(0,T;L = ) _ _^2(d-( s )-l/fc)' 

the last term is estimated the same way, for (|5.8|) still holds when replacing s with 
r: 

C 



Kccll L 5 (0,T;L") < . _ 5(r)-l/ 9 



We infer: 



thus 



.£|2\ „.e 



C 



U free\\ Lq >( T ;L>-') - , 

(l-r) 7 * 



(5.10) \\ w ^ Lx(0T . L2) <c(^^j 



7-1-i 

' Q 



Now apply the operator J e to l|5.7l) . Since J £ and U e commute, it yields, 
J £ {t)w e = U e (t)J e (0)r e - ie 1 ' 1 f W(t - s)J E (s) ((l^r 7 * |u £ | 2 ) u E ) (s)ds. 



The action of J £ on the nonlinear term is described by 1)2.9(1 . In order to apply 
Prop. 15. 31 like before, we take now 



and 



(5.11) F £ (w £ ) = (l^r 7 * |u e | 2 ) J £ {t)w £ + 2Rc(\x\-~< * {Wj e (t)w £ )y 
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The first term on the r.h.s of l|5.11[) leads to an equation which is very similar to 
(|5.7|) . with w £ replaced by J £ w £ and is treated by the same computations as above. 
For the second term, we estimate by Holder, by the Hardy-Littlewood-Sobolev 
inequality and then again by Holder: 



2Re(|zp * (Wr(t)w'))u s \\ Lri < * (WJ s (t)w s ) \\ L01 \\u s (t)\\ 



<\\u^t)\\ L 4r(t) W m)\\ Lr \\um)\\ LS 

with r, s as stated in l|5.2|) and p- = 4j- + i . Here the condition to use the Hardy- 
Littlewood-Sobolev inequality is 7 > 8(r) + S(s), which is always satisfied by Ij5.2|) . 
By applying now l|5.9|) we continue to estimate 

(| cosi|) v ; 

Then we apply, like before, Prop. 15.31 and estimate the term for h £ as above, to 
obtain: 

e 



(5.12) || J e w e \\ L ~ i0tT . <L2) < C||Vr £ || L2 + C - 



2 T 



Combining (|5.1U|) and l|5.12|l yields, along with (|2.8|) . 



Vt€[0,T], \\w e (t)\\ L „< 



c 



7-1- 



Therefore, choosing e sufficiently small and A sufficiently large, we deduce that we 
can take T = \ — Ae. With the result of Lemma lT!2l on the limit of itf rce , this yields 
Prop. 021 away from the focus, for A £ — Id and J £ . The case A £ — H £ on this 
time interval is now straightforward. 

The remaining parts of the proof for Prop. IS~2l are done as in [H] with the method 
changes as in the part shown above. It remains to show that the approximations 
in the two different regimes match at = — Ae, and that the influence of the 
harmonic potential is small near the focus so that the propagation there is given 
by 

(5.13) 1 B (t,x) = -^il>' 1 



where ip is the solution of l|1.9|) subject to the following initial condition at t = —00 

UH)^(t)| t= _ 00 = e i ^/. 

This solution exists according to Proposition ^. II 
Then the following asymptotic is proven: 



lim sup sup 

e^O f-Ae<t<f+Ae 



A e {t) (u £ (t, x)-v £ (t,x)) 



— ► 0, 

L 2 A^+oo 



with A £ (t) being one of the operators Id, J £ or H £ . Since these parts are quite 
similar to the treatment in we do not repeat them. 

After the crossing of the first focus, the solution is again propagated linearly and 
at subsequent focusing points this process is iterated. □ 
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6. Formal Computations and Discussions 

6.1. The case a = 7 = 1 (in 3-d: Schrodinger-Poisson). We saw in Sectional 

that when a = 7 > 1 , the nonlinear term in 1)1. 6J1 has a leading order influence near 
the focuses, and only in these regions. On the other hand, if a = 1 and 7 < 1, 
Section 01 shows that the Hartree term cannot be neglected away from the focuses. 
These two cases suggest that when a — 7 = 1 , the nonlinear influence is everywhere 
relevant. The aim of this final section is to give convincing arguments that this is 
the case. 

For the influence near the focuses, we need the scattering theory for (|1.9|) at 
7 = 1. In this long range scattering case, modified scattering operators are needed 
instead of the ones described in Prop. 15.11 Hayashi and Naumkin J^j obtained 
an asymptotic completeness result for n > 2 with smoothness assumptions which 
are applicable to our situation. On the other hand, they could not obtain wave 
operators. Ginibre and Velo [141 ITB] obtained modified wave operators for (|1.9|) 
with 7=1 using Gevrey spaces by a technically involved method. A drawback of 
both these results is that they include a loss in regularity. 

To show how the long range scattering theory fits into our framework we report 
(a particular case of) the result of Hayashi and Naumkin • 

Proposition 6.1 (JT6J). Assume n = 3, cp € and 5 = \\<p\\ji is sufficiently small. 
Let tp € C(R, E) be the solution of 1)1.9(1 with ipt=o = f- Then there exists a unique 
function %f> + S H a '° n H ^ , | < a < 1, such that 



— » 0, 

JJ2 t^ + OO 



\t = 



m - ex P (i (m- 1 * i^ + i 2 ) (f) log \t\) u(tty 

where H a l3 = j> e S' | ||(1+ \x\ 2 )^ 2 (l + A) a / 2 (t)\\ L 2 < 00}. 

To summarize very roughly, the results in [141 115| consist in showing that given 
some "0+ (or ip- f° r an asymptotic behavior for t — > —00), one can find ip solving 
(11.9(1 such that the above asymptotics holds. 

Analogoulsy to the treatment of long-range scattering in 0] , one can now define 
g £ (t,x) := (Ixj" 1 * |/| 2 ) (x)log(^p) (compare with ((4.11) ) and add the phase o g 
to the initial data in ((!.() (I . This yields: 

u e \ t=0 = /(^e-^M- 1 *!/! 2 )^)^. 

Using the modified scattering operators from the results of ^] we get, at 
least formally, for < t < ir/2, 

" £ (*- x ) ~ f (A) e-t^Wit,^) as e _> . 

(costy/ 2 Vcost/ 

This asymptotic also stems from the same computations as those performed in 
Section f4.1l Notice that the matching for |t — ?| = 0(e) is similar to the one 
in jSJ, except that we now have to take the presence of g e into account. This is 
where changing the integration from to t in 1)4.1)) into the above definition of g e 
makes the matching possible. Indeed, for \t — J | = 0(e), we compare u e with the 
function v e given by 1)5. 13)1 . where ip is now the solution given by the long range 
wave operators constructed in |141 115) . To make this statement more precise and 
the link between 1)4.1)1 and the definition of g £ more explicit, notice that we have, 
as t ^ ~2 

g e (t, ~ (M" 1 * l/l 2 ) (-^) 1o s (^r) ( phase shift for ^ 



2 



-rixi- 1 * 1 " 2 



f x \ f dr 

|/| 2 ) / (compare with 0~TJ) ■ 

V cost J J arccose COST 
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The effects of the nonlinearity show up in g e . Using the scaling l|5.13[) we can 
then (formally) continue with Prop. IrTTl for tt/2 < t < 3ir/2, 

I cos V / Vcost/ 

where S is the map S : V J - l— * V-'+j where ■0- is the asymptotic state of the result 
of 53 j which yields some solution ip to (jl.9(l . and t/> + is provided by Prop. 16.11 h e 
is given by 



/r 



■{t,x) := - (N- 1 * \FoSor- l f\*) (a?) log fe^i) 



The action of FoSoT^ 1 on / accounts for nonlinear effects taking place at the focus, 
and h E for nonlinear effects after the focus. So the influence of the nonlinearity will 
be relevant at all times. 

The impossibility to define a scattering operator for this case is one of the reasons 
why this argument is only formal. 

Remark. A rigorous result could be obtained with the same approach as in [3]. It 
would consist in studying the system of linear equations with a nonlinear coupling, 

ied t u £ + -e 2 Au e = ^-u £ , 
2 2 

ied t u e + \e 2 Au E = ^-u e + e (\x\~ x * |u £ | 2 ) u e . 

The first equation is solved explicitly thanks to Mehler's formula, and the second one 
is a linear Schrodinger equation with a harmonic potential and a time-dependent 
perturbation. With the oscillatory integral used in Section 0] and adapting the 
results of one could prove similar asymptotics to those stated above. 

6.2. The case of an additional local strong nonlinearity. We now consider 
equation l|1.6|l with an additional nonlinear term that is a multiplication operator 
with a power of the density |u e | 2 . 

Such equations arise in the modeling of effective one particle Schrodinger equa- 
tions where "exchange terms" like in the Hartree-Fock equation are simplified to 
functionals of the local densities, i.e. time dependent density functional theory, 
with the Schrodinger- Poisson-Xa equation as the simplest of such models (see 
and pQ for a heuristic derivation and numerical simulations). Note that the addi- 
tional "local" term has the opposite sign than the Hartree term (corresponding to 
the physical fact that the "exchange-correlation hole" weakens the direct Coulomb 
interaction). 

We will hence consider the following class of semi-classical Hartree equations 

(6.1) ied t u £ + \e 2 Au £ = ^-u e + e a * \u e \ 2 ) u e - e^\u e \ 2(T u e , 

with a > 1, j3 > 1, 7 > for x € K n , and with a a that is sub-critical with respect 
to finite time blow-up, i.e. - > a > 0. 

We can now discern the influence of the two nonlinear terms in the classical limit 
in terms of: 

• The size of the scaling exponents a, (3 with respect to the critical value . 

• The relation between the scaling and the "strength" of the nonlinearities 
determined respectively by 7 and a. 

If we take a > 1 and f3 > 1, by jH| and Section [3] we find that the classical limit 
is given by the linear propagation as long as no focusing occurs. At the focus, the 
relevant discrimination is a — (3/n or < (3/n for the power nonlinearity and 7 = a 
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or < a for the Hartree term. If a = (3/n and 7 < a, the crossing of the focus 
will be described by the scattering operator for NLS (when it is defined), if on the 
other hand a < (3/n and 7 = a (and the assumptions of Prop. [K~T1 are satisfied), 
focus crossing will be determined by the scattering operator of Prop. 15.11 If both 
nonlinearities are at the critical strength (0 — (3/n and 7 = a), then both will have 
an influence in crossing the caustic. If, on the other hand, both a < (3/n and 7 < a, 
the nonlinear influence will be negligible everywhere. 

If at least one of the scaling exponents a and j3 is equal to 1 and, at the same 
time, both a < f3/n and 7 < a, the corresponding nonlinear term will be relevant 
in the WKB propagation before the focusing. At the focus, the nonlinear terms will 
not be relevant and the crossing of the focus will be as in Prop. 14.11 If a = (3/n and 
7 = a then there will be a nonlinear influence everywhere and long range scattering 
for NLS and/or Hartree has to be taken into account. 

The influence of the nonlinear action for the single power NLS and the Hartree 
equation is summed up in two tables, for Hartree the table is given in the intro- 
duction, for single power nonlinear Schrodinger equation, it is stated in 3 . The 
behavior of l|6.1|l can be described by independently superposing these two tables. 
The following table is an extract from that superposition: 





a > 7 and (3 > an 


a = 7 or (3 — an 


a > 1 and (3 > 1 


Linear WKB, 


Linear WKB, 




linear focus 


nonlinear focus 


a = 1 or (3 = 1 


Nonlinear WKB, 


Nonlinear WKB, 




linear focus 


nonlinear focus 



"Nonlinear WKB" respectively "nonlinear focus" here stands for an influence from 
at least one of the nonlinear terms away from the focus or close to the focus. 

6.3. Wigner measures. We already mentioned in the introduction the work of 
Zhang, Zheng and Mauser |31| where the (semi) classical limit of the Schrodinger- 
Poisson equation with no smallness assumption (on the initial data or the nonlin- 
earity) is studied by means of Wigner measures. Wigner measures have proven to 
be efficient tools for linear semi-classical problems and for homogenization limits; 
see |26| for an overview on Wigner measure limits of Hartree equations. Wigner 
measures have the merit that in phase space the caustics of physical space are some- 
what unfolded and that generally, results globally in time are possible. 
In |S] , the Wigner measure of the nonlinear Schrodinger equation with power-like 
nonlinearity studied in [3] is investigated. It is shown that the Wigner measure 
leads to an ill-posed problem whenever nonlinear effects at the focal points come 
into play. In other words, the Wigner measure can only be valid as long as no 
caustic appears. We briefly discuss the Wigner measures of 1)1.6(1 in view of these 
results. 

The Wigner measure of the family (u £ (i))o<e<i, which is bounded in L 2 , is the 
weak limit under e — > (up to an extraction) of its Wigner transform, 

This limit is a positive radon measure /1 and is in general not a unique limit. 

- linear case: Case a > 7, a > 1: 
By the result of Prop. 13. II and the asymptotics of U[ ICC in Lemma f4. 21 the Wigner 
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measure 



[a for t < 7r/2 of the family (ii E (i))o<e<i is 



2 

dx ® (5. 



x tan i ■ 



For ir/2 < < < 7r, the Wigner measure of {u £ (t))o <e <i (denoted by fj, + ) is the same: 
/i + (i,x,^) = n~(t,x,£). At t = n/2, the limits from above and below are: 
fi~(t,x,g) =Um t ^ 7r/2 +/i+(t,x,0 = |/(£)| 2 d£ ® 

- nonlinear WKB, linear focus: Case 7 < a = 1. 

The asymptotics of u £ are stated in Prop. 14.11 The additional phase term g is of 
order 1 and does not change the Wigner measure of (u e (t))o <e <i, so in this case fx" 
and n + are the same as in the previous case: the Wigner measure does not "see" 
the nonlinear effect g . 

- linear WKB, nonlinear focus: Case 7 = a > 1. 

The asymptotics of Prop. l5.2l involve. for t > ir/2, the scattering operator associated 
with the unsealed equation (|1.9(l . For t < 7r/2, the Wigner measure of (it E (i))o<e<i 
is still the same as above, but for ir/2 < t < ir, we have 

2 



^ v J cost ™ V 



dx (8) Extant, 



v COS t / 

where S 1 is the scattering operator for l|1.9|) and the Fourier transform. 

- nonlinear WKB, nonlinear focus: Case 7 = a = 1. 
The asymptotics for this case of (the formal computation) Prop. 16.11 include an 
additional phase term which is of order log e and a modification of the initial data 
of the same order of magnitude. Both do not alter the Wigner measure, since 
they are dominated by the scaling of the Wigner transform, and thus the Wigner 
measure is the same as in the previous case. 
For the last two cases, the limits at t = tt/2 are 

lim /i~(*,*,o = i/(Oi a de 

f— >7i72- 

6-2 ' , , 

lim ii + (t,x,£) = \FoSof- 1 f(n\ 2 dt<g>5(x). 

t^7r/2+ 

The idea of 5 is to find now two profiles f\ and fa for which |/i| 2 = I/2I 2 , but at 
the same time \T o S o T^ 1 f\\ 2 ^ | T o S o J-~ 1 f-2\ 2 ■ Then the Wigner measures 
of the corresponding families (u E (i))o< £ <i, j — 1,2, will be equal up to the focus, 
but different after the focus, i.e. /i^ = [1% but /xj 7^ ^ . So after the caustic point 
the Wigner measure will not be unique anymore in the case where the nonlinearity 
is relevant at the focus. These profiles were constructed using an expansion of S 
around the origin. Since our problem is very similar to the one studied there, we 
expect a similar result to hold for equation i|1.6fl . he. we expect the Wigner measure 
to lead to an ill-posed problem if there is a nonlinear influence at the caustic. 

In view of the result of |3Xj , note that the non- uniqueness of the weak solutions for 
Vlasov-Poisson with measures as initial data and the non-uniqueness of the Wigner 
measure of a given e-dependent family of solutions coincide, such that there is no 
contradiction with the global and unique semi-classical limits of the Hartree type 
equations obtained here. 
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